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Abstract 


In  a  comparison  of  measured,  simulated,  and  theoretical  calculations  of 
radar  cross  section  (RCS),  the  radar  target  was  a  simple  missile  shape  (a 
cylinder  with  hemispherical  ends).  For  the  measurements,  the  target  was 
rotated  on  a  turntable  from  -5°  to  365°  in  azimuth,  and  its  RCS  was  measured 
with  a  W-band  frequency-agile  instrumentation  radar.  (Clutter  and  noise  in 
the  measured  RCS  data  were  reduced  by  spatial  fQtering.)  For  modeling  the 
target  geometry,  the  ACAD  (advanced  computer-aided  design)  geometric 
modeling  program  was  used  to  create  both  an  IGES  (initial  graphics  ex¬ 
change  specification)  and  a  facet  model.  For  modeling  the  RCS,  three 
high-frequency  prediction  codes  were  used  (Xpatch,  CADDSCAT,  and 
NcPTD),  as  well  as  a  point  scatter  model.  Various  methods  were  used  to 
process  and  simulate  the  target  returns,  depending  upon  the  aspect  angle. 
Agreement  between  the  measured  and  simulated  bandwidth-averaged  RCS 
values  depended  on  the  portion  of  the  target  that  dominated  the  radar 
return:  agreement  was  close  for  radar  returns  dominated  by  returns  from  the 
hemispherical  ends  of  the  target.  However,  for  aspect  angles  near  broadside 
to  the  target  (for  which  the  cylindrical  part  of  the  target  dominated  the 
returns),  RCS  measurements  suggested  an  interaction  between  the  target 
and  the  pylons  supporting  it. 
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1.  Introduction 


The  proliferation  of  tactical  ballistic  missiles  throughout  the  world  has  pro¬ 
vided  motivation  for  the  United  States  to  develop  theater  ballistic  missile 
defense  systems  such  as  Patriot,  Theater  High-Altitude  Air  Defense 
(THAAD),  Navy  Area,  Navy  Theater  Wide,  and  Medium-Extended  Air 
Defense  System  (MEADS).  Previously  developed  missiles  using  conven¬ 
tional  fuzing  systems  have  difficulty  destroying  fast-moving  "Scud-like" 
missiles.  To  defeat  targets  with  high  speed  and/ or  agility  and  intensive 
countermeasures,  new  missiles  are  being  developed  that  use  techniques 
such  as  hit-to-kill  or  guidance  integrated  fuzing  (GIF).  Because  of  the  high 
cost  of  field  testing,  computer  simulations  will  have  a  significant  role  in  the 
development  of  these  systems. 

Important  to  such  computer  simulations  are  validated  and  com¬ 
putationally  efficient  radar  backscatter  target  models.  There  are  several 
validated  radar  cross  section  (RCS)  predictive  codes,  but  their  levels  of 
validation  depend  upon  target,  frequency,  range,  and  computer  processing 
time  [1].  Several  codes  perform  well  when  used  for  a  specific  application, 
but  no  code  has  been  found  that  is  computationally  efficient  and  can 
model  all  the  features  desired  for  simulating  GIF  for  complex  targets  (these 
include  multi-bounce  effects,  near-field  effects,  variable  antenna  patterns, 
multiple  targets,  and  traveling  waves).  This  report  compares  measured, 
simulated,  and  theoretical  calculations  of  bandwidth-averaged  RCS  for  a 
missile-shaped  target  at  W-band. 

The  Army  Research  Laboratory  (ARL)  measured  the  radar  backscatter  of  a 
test  target  consisting  of  a  cylinder  capped  by  two  hemispheres,  as  the  tar¬ 
get  rotated  on  a  turntable.  Measurements  were  made  at  Aberdeen  Proving 
Ground  with  a  frequency-agile,  coherent,  polarimetric,  instrumentation  ra¬ 
dar  built  by  ARL.  The  measurements  were  made  in  the  far  field  of  the  tar¬ 
get,  except  for  aspect  angles  near  broadside  to  the  cylindrical  portion  of  the 
target.  Although  the  full  linear  polarization  scattering  matrix  was  meas¬ 
ured,  this  report  analyzes  only  the  vertical-transmit,  vertical-receive  (W) 
polarization. 

The  radar  returns  were  simulated  with  three  mature  high-frequency  RCS 
prediction  codes*  (Xpatch  [2],  CADDSCAT  [3],  and  NcPDT  [4])  and  a  point 
scatter  model  (PSM).  These  codes  use  physical  optics  (PO)  and  the  physical 
theory  of  diffraction  to  calculate  electric  fields  and  RCS  [5].  They  caimot 
simulate  traveling  or  creeping  waves. 

Techniques  were  developed  to  reduce  and  separate  the  noise  and  clutter  from 
the  target  signature.  To  enable  a  direct  comparison  between  measured  data 
and  simulated  results,  I  processed  the  measured  data  with  spatial  filters  that 
reduced  clutter,  noise,  and  (if  present)  traveling  and  creeping  waves.  For  most 
aspect  angles,  the  scattering  due  to  traveling  and  creeping  waves  and  noise 
appeared  to  be  negligible.  The  RCS  at  aspect  angles  where  such  scattering 
may  have  been  significant  was  calculated,  but  not  analyzed. 


*Xpatch  and  NcPDT  are  a  product  ofDEMACO,  and  CADDSCAT  is  a  product  ofBoeing/McDonnell-Douglas.  The 
federal  government  has  contributed  to  the  development  of  these  codes  through  the  Electromagnetic  Code  Consortium. 
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2.  Modeling  and  Experimental  Setup 

2.1  Target  Characteristics 

The  test  target  was  a  cylinder  with  hemispherical  ends,  as  pictured  in  fig¬ 
ure  1.  The  cylinder  is  made  up  of  three  sections,  which  were  threaded  at 
the  ends  so  that  they  could  be  screwed  together.  The  cylinder  is  180  in. 
long,  with  each  section  being  60  in.  long.  The  radius  of  the  hemispheres  is 
4  in.,  so  that  the  diameter  of  the  cylinder  is  8  in.  The  hemispheres  are  also 
threaded  so  that  they  can  be  screwed  onto  the  ends  of  the  cylinder,  and 
there  is  a  3/8-in.  extension  on  each  hemisphere  for  the  threads.  The  gaps 
between  the  seams  of  the  cylinder  were  all  less  than  0.002  in.  The  overall 
length  of  the  test  target  is  188.75  in.:  180  in.  for  the  cylinder,  4  in.  for  each 
hemisphere,  and  0.375  in.  for  a  cylindrical  extension  on  each  hemisphere. 
The  approximate  weight  of  the  target  is  400  lb.  (Fig.  3,  sect.  2,  illustrates  the 
dimensions  of  the  target.) 

The  target  was  supported  by  two  buoyant  Styrofoam  pylons  that  were 
shaped  to  have  a  minimum  RCS  at  aspect  angles  near  0°  azimuth.  The 
width  of  the  top  of  the  pylons  (shown  in  fig.  1)  was  approximately  1/8  to 
1/4  of  an  inch,  and  increased  to  approximately  1.5  in.  at  the  bottom  of  the 
test  target.  These  widths  are  important,  because  a  radar  signal  would  be 
slightly  delayed  while  propagating  through  the  pylons  towards  a  target 
for  certain  aspect  angles.  The  exact  dimensions  of  the  pylons  were  difficult 
to  measure  and  quantify  because  they  are  curved  three-dimensional  (3-D) 
objects. 

During  the  measurements,  the  target  was  placed  30  in.  from  the  center  of 
rotation  of  the  turntable,  with  0°  azimuth  angle,  corresponding  to  the  ge¬ 
ometry  shown  in  figure  2.  The  depression  angle  was  1.45°  with  respect  to 
the  radar  (a  positive  depression  angle  implies  that  the  radar  was  above  the 
target).  The  peak  of  the  antenna  pattern  was  pointed  towards  the  target  at 

Figure  1.  Test  target 
and  pylons  on 
turntable  (aspect  angle 
near  broadside). 
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Figure  2.  Geometry  of 
test  target  at  0° 
azimuth  angle. 
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0°  azimuth  angle.  The  two-way,  3-dB  width  of  the  antenna  pattern  was  ap¬ 
proximately  4.4°,  which  was  slightly  smaller  than  the  target  at  broadside  at 
a  range  of  46.6  m.  The  exact  location  of  the  target  and  the  antenna  pattern 
were  determined  through  a  Monte  Carlo  simulation  based  upon  the  radar 
measurements. 


2.2  Target  Model 

For  input  to  the  NcPDT,  CADDSCAT,  and  Xpatch  simulation  codes,  I  gen¬ 
erated  the  test  target  surface  using  the  ACAD  (advanced  computer-aided 
design)  software  program.  All  the  simulation  codes  could  handle  multiple 
geometric  representations;  I  chose  two  commonly  used  representations  for 
the  target  model;  an  ICES  (initial  graphics  exchange  specification)  format 
and  a  facet  model  format.  I  created  arbitrary  surface  patches  by  revolving 
lines  and  curves;  to  determine  the  number  of  surface  patches  to  use,  I  in¬ 
creased  the  number  of  patches  until  the  RCS  converged  to  a  constant  for 
the  hemisphere  portion  of  the  target.  (This  approach  is  based  on  the  fact 
that  the  RCS  of  a  sphere,  in  the  far  field,  is  constant  and  independent  of  as¬ 
pect  angle  at  high  frequencies.)  I  calculated  the  standard  deviation  of 
bandwidth-averaged  RCS  simulated  with  CADDSCAT  at  different  aspect 
angles  and  compared  these  values  to  a  threshold.  I  then  increased  the 
number  of  surface  patches  in  the  IGES  representation  imtil  the  standard 
deviation  of  the  simulated  RCS  was  less  than  or  equal  to  0.1  dB.  This  re¬ 
sulted  in  a  model  for  the  hemispheres  that  contained  18  spherical  slices  of 
20°. 

I  followed  a  similar  procedure  for  Xpatch,  using  the  same  ACAD  model, 
but  with  the  facet  model  representation.  For  CADDSCAT,  I  implemented  a 
rectangular  gap  model  using  spline  representations  of  the  gaps  created  by 
the  physically  segmented  target. 

The  cylinder  section  of  the  target  was  modeled  with  six  smaller  cylinders 
as  shown  in  figure  3.  The  lengths  of  the  cylinders  were  determined  by  the 
starting  and  ending  horizontal  position  of  the  top  of  the  pylons.  The  four 
larger  cylinders  each  contained  18  cylindrical  slices  of  20°.  The  two  smaller 
cylinders  contained  16  cylindrical  slices  of  20°  and  4  of  10°.  These  smaller 
cylinders  were  divided  into  top  and  bottom  sections.  (In  the  figure,  the 
bottom  cylinders  are  offset  from  the  top  cylinders  for  visual  clarity.)  The 
length  of  the  bottom  sections  corresponded  to  the  measured  length  of  the 
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Figure  3.  Geometric 
target  model  8  in. 

represented  by  IGES 
128  surface  patches. 

Smaller  half  cylinders 
at  bottom  are  offset 
from  main  target  for 
clarity. 

Styrofoam  columns  supporting  the  target  (and  obscuring  parts  of  it);  the 
height  of  the  half  cylinder  was  set  to  4  in.,  the  position  of  specular  reflec¬ 
tion  for  a  radar  at  a  0°  depression  angle.  However,  the  measurements  indi¬ 
cated  that  on  the  actual  target,  this  height  varied  slightly  from  4  in.  Since 
the  simulated  RCS  was  very  sensitive  to  the  pylon  vertical  location  for  val¬ 
ues  near  4  in.,  the  inability  to  input  the  true  dimension  affected  the  model 
results.  Figure  3  shows  the  ACAD  model  of  the  test  target. 

For  diagnostic  purposes,  a  PSM  was  constructed  to  simulate  the  returns 
from  the  hemispheres.  The  model  consisted  of  a  single  scatterer  with  a  con¬ 
stant  amplitude,  with  a  position  that  tracked  the  specular  reflection  point 
for  the  appropriate  hemisphere. 

2.3  Radar  Characteristics 

ARL  has  a  coherent  fully  polarimetric  radar  system  capable  of  operation  at 
X-band,  X^-band,  or  W-band  [6].  Table  1  lists  the  current  characteristics  of 
the  W-band  radar  system  that  was  used  to  measure  the  test  target. 

The  in-phase  and  quadrature  (I&Q)  outputs  of  the  system  are  digitized  in 
the  radar  to  12-bit  accuracy  for  each  pulse.  Four  analog  to  digital  (A/D) 
converters  sample  the  vertical  and  horizontal  I&Q  returns  from  eight  verti¬ 
cally  transmitted  pulses,  followed  by  the  returns  from  eight  horizontally 
transmitted  pulses  for  each  of  the  256  steps  in  frequency.  Data  collection  is 
triggered  when  a  preselected  angle  of  rotation  has  been  reached.  The  data 
are  then  demultiplexed,  averaged  over  eight  pulse  groups,  and  stored  as 
256  I&Q  pairs  for  both  the  vertical  and  horizontal  receive  channels.  Time 
and  turntable  angle  are  also  recorded  with  each  ramp  of  data. 

2.4  Measured  Data 

The  target  measurements  were  made  on  9  June  1997,  at  the  Aberdeen  Test 
Center  (ATC)  outdoor  high-resolution  millimeter-wave  (MMW)  signature 
measurement  facility  at  Airbase  Range  8.  The  target  was  rotated  370°;  this 
range  includes  an  extra  5°  on  either  side  of  the  0°  mark  at  the  start  and  end 
of  the  rotation,  which  was  to  be  used  for  checking  signal  level  drift  during 
the  measurements.  The  raw  data  were  calibrated  with  software  developed 
by  ARL  [6].  The  data  consisted  of  I&Q  values  for  256  frequencies  incremen¬ 
tally  stepped  between  93.488  and  95.0  GHz  for  four  polarizations;  these 
I&Q  values  were  sampled  every  0.006°.  A  ramp  of  data  is  defined  as  the 
I&Q  values  associated  with  all  the  linearly  swept  frequencies  for  a  single 
polarization  and  aspect  angle.  The  calibrated  data  were  coherently 
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Table  1.  Radar 
characteristics  at  W- 
band. 


Parameter 

Value 

Maximum  jfrequency 

95.0  GHz 

Minimum  frequency 

93.488  GHz 

Rf  bandwi(Jth 

1511.64  MHz 

Frequency  step 

5.928  MHz 

Peak  transmitted  power 

+13  dBm 

Pulse  width 

100  ns 

Pulse  repetition  frequency 

1.0  MHz 

Transmitted  polarization 

VorH 

Received  polarization 

VandH 

3-dB  beam  width  (two  way) 

4.4° 

System  noise  figure 

11  dB  single  sideband 

Minimum  detectable  signal 

-62  dBsm 

Polarization  isolation 

30  dB 

averaged  over  1024  ramps,  or  approximately  6°  segments,  to  form  an  I&Q 
offset.  This  offset  was  subtracted  from  each  ramp  of  data  in  the  segment. 


To  ensure  that  the  95-GHz  radar  was  operating  properly,  I  performed  an 
error  analysis  on  the  calibration  measurements,  determining  the  ampli¬ 
tude  and  phase  stability  of  a  trihedral  and  a  dihedral  on  a  rotating  mount. 
Figures  4  and  5  show  the  RCS  and  the  variation  in  phase  as  a  function  of 
frequency  for  co-polarization  for  the  dihedral  rotated  through  22.5°.  I  cal¬ 
culated  the  variation  in  phase  by  subtracting  the  phase  associated  with  the 
frequency  with  index  j  from  the  phase  of  frequency  with  index  j  +  1,  then 
properly  unwrapping  the  phase  difference.  I  then  chose  the  22.5°  dihedral 
for  analysis  because  it  has  a  minimal  impact  on  calibration.  After  each  fre¬ 
quency  was  averaged  over  2000  measurements,  I  found  that  the  high  fre¬ 
quencies  have  a  lower  signal  to  noise  ratio  (SNR).  Further  analysis  indi¬ 
cated  that  only  the  last  20  frequencies  were  of  poorer  quality;  these  were 
eliminated  from  further  analysis.  Table  2  shows  the  standard  deviation  of 
the  RCS  and  phase  calculated  over  frequency  with  the  last  20  frequencies 
eliminated.  The  small  values  indicate  that  the  radar  is  functioning  properly 
over  the  reduced  bandwidth.  Table  2  provides  a  lower  boimd  for  the  meas¬ 
urement  error. 


2.5  Simulation  Models 

The  radar  returns  were  simulated  with  Xpatch,  CADE>SCAT,  NcPDT,  and  a 
PSM.  The  simulations  were  run  with  version  2.1  of  Xpatch  (with 
frequency-domain  analysis),  beta  version  1.6.6  of  CADDSCAT,  and  version 
1.2  of  NcPDT  with  the  program  Ncl28.f  (Ncl28.f  required  some  small 
modifications).  The  geometric  models  input  to  the  simulations  were  gener¬ 
ated  with  ACAD  (version  9.0c),  a  geometric  modeling  program  developed 
by  Lockheed  Martin.  I  used  the  PSM  to  validate  the  processing  software 
and  to  estimate  the  processing  errors. 

Each  code  has  different  capabilities.  CADDSCAT  and  NcPDT  simulations 
can  be  performed  in  both  the  near  and  far  fields,  while  Xpatch  simulations 
can  be  performed  only  in  the  far  field.  In  NcPDT,  simulations  can  be  per¬ 
formed  with  a  measured  or  simulated  antenna  pattern.  Both  CADDSCAT 
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Figure  4.  RCS  for  a 
22.5°  dihedral 
averaged  over  2000 
measurements. 
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Figure  5.  Phase 
variation  for  a  22.5° 
dihedral  averaged 
over  2000 
measurements. 
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and  Xpatch  allow  multiple-bounce  simulations.  None  of  the  codes,  how¬ 
ever,  can  simulate  traveling  or  creeping  waves. 

Par-field  simulated  data  were  generated  with  Xpatch,  CADDSCAT,  and  a 
PSM  that  attempted  to  duplicate  the  measured  data.  The  simulated  aspect 
angles  and  frequencies  matched  the  measurement  parameters.  However, 
the  simulated  antenna  pattern  was  uniform,  and  the  range  was  infinite. 
Because  of  symmetry,  the  simulations  were  run  only  from  0°  to  90°  in  azi¬ 
muth  angle.  I  ran  both  predictive  codes  using  a  single-bounce  calculation 
for  a  perfect  electric  conductor  (PEC)  with  no  edge  diffraction  or  shadow¬ 
ing.  The  default  integration  densities  were  used.  CADDSCAT  was  initially 
run  with  gap  diffraction  included,  but  the  maximum  contribution  to  the 
RCS  was  below  -40  dBsm,  so  gap  diffraction  was  subsequently  omitted. 
Initially,  a  facet  divergence  file  was  included  in  the  Xpatch  simulation,  but 
since  it  did  not  significantly  affect  the  RCS,  it  was  subsequently  omitted. 

Near-field  simulated  data  were  generated  with  CADDSCAT  and  NcPDT, 
and  input  into  an  RCS  model  that  included  a  target/ pylon  interaction.  The 
simulated  range  corresponded  with  the  measurement  range,  but  I  reduced 
the  number  of  frequencies  and  aspect  angles  to  limit  computation  time.  A 
single  frequency  was  simulated  for  aspect  angles  near  90°  azimuth  angle. 
In  CADDSCAT,  I  used  a  uniform  antenna  pattern,  and  in  NcPDT,  I  used 
the  measured  antenna  pattern.  The  PO  calculations  in  NcPDT  were  per¬ 
formed  with  the  first-order  method  and  the  default  integration  density, 
and  the  antenna  pattern  was  assumed  to  be  roU  S5mimetric  [3]. 

For  CADDSCAT  and  NcPDT,  the  target  model  was  represented  in  ICES 
128  format,  and  for  Xpatch  in  DEMACO  facet  model  format.  For 
CADDSCAT  and  NcPDT,  each  hemisphere  and  each  complete  cylinder 
contained  18  patches.  For  Xpatch,  the  target  model  contained  10,440  facets, 
and  each  hemisphere  contained  4561  facets.  The  facet  model  was  gener¬ 
ated  from  the  IGES  128  model  by  the  conversion  utilities  bundled  with 
Xpatch.  The  tolerance  selected  to  generate  the  facets  was  0.003  in.,  or  0.074 
mm;  this  dimension  is  much  smaller  than  the  average  wavelength  used  in 
the  simulation  (2.3  mm). 


Table  2.  Standard 
deviation  of  RCS  and 
phase  variation  of 
22.5°  dihedral 
calculated  over 
frequencies  that  have 
a  high  SNR. 
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3.  Procedure 


3.1 


Processing 


I  processed  the  radar  data  to  reduce  clutter  and  noise,  correct  for  the  effects 
of  the  antenna,  and  fine  tune  the  RCS  range  correction.  Figure  6  shows  a 
block  diagram  of  the  processing  procedure.  Variables  with  capital  letters 
represent  signals  in  the  frequency  domain,  and  variables  with  small  letters 
represent  signals  in  the  spatial  domain.  The  input  signal  S  was  loaded, 
consisting  of  I&Q  data  for  different  frequencies  and  aspect  angles.  A  Kaiser 
window  (IV)  then  filters  the  signal  to  reduce  blurring  in  the  subsequent 
processed  image.  Next,  a  discrete  Fourier  transform  (DFT)  was  performed 
to  image  the  data,  and  then  a  bandpass  (B)  filter  was  applied  to  remove 
clutter  and  noise  outside  the  target  region.  An  inverse  was  performed 
to  provide  I&Q  values,  and  an  inverse  of  the  original  Kaiser  window  (W) 
was  applied.  The  transfer  function  is  given  by 


H{z)  = 


(SW)  *  b 
SW  ' 


(1) 


where  S  is  the  I&Q  data,  W  is  the  coefficients  for  the  Kaiser  window,  b  is 
the  inverse  DFT  of  the  bandpass  filter  B,  and  the  operator  denotes 
convolution. 


The  procedures  used  to  process  the  data  were  get,  filter,  image,  extract,  in¬ 
verse,  and  bandwidth  average,  as  shown  on  the  left  side  of  the  flow  chart  in 
figure  7.  The  procedures  were  implemented  as  hmctions  in  MATLAB  (ver¬ 
sion  4.2).  The  major  components  of  the  procedures  are  shown  on  the  right 
side  of  the  figure.  The  components  are  divided  between  1-D  processing 
techniques  (which  I  refer  to  collectively  as  the  high-range  resolution  (HRR) 
processing  technique)  and  2-D  processing  techniques  (which  I  refer  to  col¬ 
lectively  as  the  inverse  synthetic  aperture  radar  (ISAR)  processing  tech¬ 
nique).  Several  components  shown  in  figure  7  are  used  by  both  techniques. 


3.1.1  GetData 

In  the  "get  data"  procedure,  l&Q  data  were  read  from  measured  and  simu¬ 
lated  data  files,  or  generated  by  a  PSM.  This  step  was  simplified  by  the  use 
of  a  common  binary  format  for  data  storage.  PSM  data  were  generated 
separately. 


3.1.2  Filter 

In  the  filter  procedure,  the  data  were  filtered  by  a  Kaiser  window  with  P  = 
1,  where  the  P  parameter  affects  the  sidelobe  attenuation.  If  the  HRR 
processing  technique  was  selected,  the  Kaiser  window  was  1-D,  and  if  the 
ISAR  processing  technique  was  selected,  the  Kaiser  window  was  2-D. 
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Figure  6.  Signal  flow  diagram  of  data  processing. 
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Figure  7.  Flow  chart  for  computing  bandwidth-averaged  RCS:  (left)  procedures  and  (right)  major 
components. 


3.1.3  Image 

In  the  image  procedure,  a  1-D  fast  Fourier  transform  (FFT)  was  performed 
to  generate  an  HRR  profile,  or  a  2-D  FFT  was  performed  to  generate  an 
ISAR  image. 

3.1.4  Extract 

In  the  extract  procedure,  the  algorithm  extracted  the  target  from  clutter 
and  noise  by  clipping  regions  of  the  HRR  profile  or  ISAR  image.  Initially, 
the  extraction  region  was  determined  by  target  geometry  and  aspect  angle. 
The  image  (or  equivalently,  the  data)  was  divided  into  two  regions,  near 
broadside  and  head  on,  for  which  different  extraction  algorithms  were  em¬ 
ployed.  The  near-broadside  region  was  defined  as  within  ±15°  of  the 
broadside  of  the  cylinder  portion  of  the  target,  and  the  head-on  region  was 
defined  as  within  ±75°  of  the  apex  of  either  hemisphere.  In  the  head-on  re¬ 
gion,  the  RCS  was  dominated  by  returns  from  the  facing  hemisphere,  and 
in  the  near-broadside  region,  it  was  dominated  by  returns  from  the  cylin¬ 
der  portion  of  the  test  target.  In  practice  the  division  between  the  two  re¬ 
gions  is  not  smooth,  so  there  will  be  a  small  discontinuity  in  the  processed 
RCS. 

The  extraction  algorithm  used  varied  with  the  processing  technique  (ISAR 
or  HRR)  and  image  region.  When  the  ISAR  processing  technique  was  ap¬ 
plied  in  the  head-on  region,  the  extraction  algorithm  was  based  upon  a 
diamond-shaped  area  that  was  translated  and  scaled.  The  center  of  the  dia¬ 
mond  corresponded  to  the  location  of  the  specular  reflection  from  the  ap¬ 
propriate  hemisphere  of  the  target.  The  size  of  the  diamond  was  empiri¬ 
cally  determined,  and  it  was  scaled  as  a  function  of  aspect  angle.  The  shape 
of  the  diamond  was  not  corrected  for  curvature,  but  the  location  of  the  cen¬ 
ter  of  the  diamond  with  respect  to  the  ISAR  image  was  corrected. 

For  ISAR  processing  in  the  near-broadside  region,  the  extraction  algorithm 
was  based  upon  a  rectangle  that  was  translated,  rotated,  scaled,  and  cor¬ 
rected  for  curvature.  The  rectangle  corresponded  to  the  geometry  of  the 
target,  but  was  expanded  and  shifted  downrange.  An  additional  rectangu¬ 
lar  region  was  included  for  aspect  angles  that  had  a  large  specular  return 
from  the  cylinder  portion  of  the  target.  To  encapsulate  the  sidelobes  in  the 
image,  this  area  was  added  for  aspect  angles  ±2.5°  from  broadside.  The 
center  of  the  rectangle  corresponded  to  the  location  of  the  specular  return 
on  the  target. 

For  HRR  processing  in  both  regions,  the  extraction  algorithm  was  based 
upon  the  target  occupying  a  series  of  HRR  cells.  The  center  of  the  extrac¬ 
tion  region  was  calculated  from  the  geometry  of  the  target.  The  extraction 
procedure  was  modeled  as  a  bandpass  filter  in  the  spatial  domain,  as 
shown  in  figure  7. 

Figures  8  and  9  illustrate  the  different  extraction  regions  for  ISAR  images 
at  azimuth  angles  of  0°  and  88°  and  an  elevation  angle  of  1.4°,  for  W  po- 


Figure  8.  ISAR  image 
of  target  for  region 
dominated  by  a 
hemisphere; 
diamond-shaped 
outline  encloses  target 
RCS. 


Figure  9.  ISAR  image 
of  target  for  region 
dominated  by 
cylinder;  trapezoid- 
shaped  outline 
encloses  target  RCS. 


larization.  The  dotted  lines  in  these  figures  denote  the  extraction  region. 
The  downrange  and  crossrange  resolution  is  approximately  0.1  m.  The 
contour  line  levels  range  from  10  to  -40  dBsm,  in  10-dB  intervals.  The  most 
prevalent  contour  level  is  -40  dBsm. 

Occasionally,  significant  blurring  occurred  in  the  image,  which  required 
the  extraction  region  to  be  adjusted.  The  most  likely  cause  was  wind 
shaking  the  target,  but  meteorological  data  were  not  available  to  confirm 
this  hypothesis.  If  the  ISAR  processing  technique  was  used  with  the  target 
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in  the  head-on  region,  and  the  computed  standard  deviation  was  greater 
than  2.4  dB,  then  the  size  of  the  extraction  region  was  increased  and  the 
processing  procedure  was  repeated.  The  size  of  extraction  region  was  in¬ 
creased  only  once  per  ISAR  image.  I  chose  the  threshold  for  the  standard 
deviation  by  examining  the  processed  data. 

Next,  if  appropriate,  I  applied  an  antenna  pattern  correction  to  the  data. 
The  method  used  depended  on  the  processing  technique  and  region.  If 
ISAR  processing  was  used  and  the  target  was  in  the  near-broadside  region, 
then  each  cell  in  the  image  was  divided  by  the  antenna  pattern  value  corre¬ 
sponding  to  its  angular  position.  If  the  target  was  m  the  head-on  region,  all 
the  extracted  cells  in  the  image  were  divided  by  a  single  antenna  pattern 
value  corresponding  to  the  geometric  location  of  the  reflection  point  from 
the  appropriate  hemisphere. 

For  HRR  processing,  the  same  antenna  pattern  correction  technique  was 
used  for  the  head-on  region.  But  if  the  target  was  in  the  near-broadside  re¬ 
gion,  the  target  was  distributed  in  crossrange,  so  the  effect  of  the  antenna 
could  not  be  removed  in  this  way. 

The  last  step  in  the  extract  procedure  was  a  fine  range  correction  of  the 
RCS.  Initially,  the  calibrated  data  were  range  corrected  to  the  center  of  the 
target,  46.66  m.  As  the  target  aspect  angle  changed,  the  range  to  the  specu¬ 
lar  reflection  point  on  the  hemispheres  changed.  This  change  was  signifi¬ 
cant,  since  it  could  cause  the  RCS  to  change  by  as  much  as  0.8  dB. 

3.1.5  Inverse 

In  the  inverse  procedure,  the  modified  image  was  converted  back  to  I&Q 
values  by  a  1-D  or  2-D  inverse  FFT  algorithm.  I  removed  the  Kaiser  win¬ 
dow  by  multiplying  the  I&Q  values  by  the  inverse  of  the  filter  coefficients, 
then  normalizing  the  RCS  so  that  power  was  conserved. 


3.1.6  Bandwidth  Average 

In  the  bandwidth-averaging  procedure,  the  I&Q  values  were  converted  to 
RCS,  then  averaged  over  frequency.  If  the  ISAR  technique  was  used  to 
process  the  data,  bandwidth-averaged  RCS  values  were  truncated  at  the 
beginning  and  ending  aspect  angles  of  the  ramps  of  data  that  were  used  to 
generate  an  image.  This  step  was  necessary  to  minimize  the  effect  of  the 
convolution  performed  in  equation  (1).  I  calculated  the  number  of  values 
truncated  by  taking  the  number  of  ramps  used  to  generate  the  ISAR  image, 
then  dividing  by  8.  For  example,  if  32  ramps  were  used  to  generate  the 
ISAR  image,  then  the  first  and  last  four  RCS  values  were  eliminated.  The 
RCS  values  at  the  remaining  aspect  angles  were  normalized  so  that  the 
same  average  RCS  was  maintained. 


3.2  Antenna  Pattern 

A 1-D  measured  antenna  pattern  was  calculated  from  a  raster  scan  of  a  re¬ 
flector  at  a  range  of  47  m.  First  the  data  were  filtered  along  the  0°  elevation 
scan  by  a  2-D  finite-duration  impulse  response  (FIR)  smoother.  Then,  the 
data  were  clipped  and  fit  to  a  second-degree  pol)momial.  The  filtered  an¬ 
tenna  pattern  was  normalized  by  its  maximum  value,  and  its  angle  look¬ 
up  table  was  offset  to  the  peak  location  of  the  fitted  polynomial.  Figure  10 
shows  the  resulting  two-way,  1-D  antenna  pattern  for  VV  polarization.  The 
two-way  beam  width  is  approximately  4.4°.  This  is  smaller  than  the  angu¬ 
lar  extent  of  the  cylinder  portion  of  the  target  at  broadside,  which  is  5.7° 
The  peak  of  the  antenna  pattern  was  offset  from  0°  azimuth  angle  by  0.66°. 

3.3  Parameters 

I  evaluated  the  effects  of  varying  several  parameters  used  in  the  HRR  and 
ISAR  processing  techniques,  performing  the  analysis  over  a  subsection  of 
the  head-on  region  to  increase  accuracy.  (I  determined  the  region  by  exam¬ 
ining  the  processed  RCS  results.)  Figure  11  shows  a  graph  of  the  unproc¬ 
essed  and  processed  bandwidth-averaged  RCS  obtained  by  the  128-ramp 
ISAR  processing  technique.  The  plot  of  the  improcessed  RCS  indicates  that 
the  noise/clutter  in  the  head-on  region  is  approximately  -3  dBsm,  which  is 
significantly  higher  than  the  predicted  RCS  of  approximately  -15  dBsm. 
The  low-noise,  head-on  region  was  determined  to  be  between  azimuth 
angles  of  0°  to  60°,  120°  to  240°,  and  300°  to  360°,  as  illustrated  by  the  verti¬ 
cal  lines  in  figure  11.  For  aspect  angles  near  the  transition  between  the 
head-on  and  the  near-broadside  region,  an  increase  in  the  RCS  was 
observed.  The  discontinuity  between  the  two  regions  was  caused  by  an  in¬ 
crease  in  the  noise  (due  to  the  use  of  a  larger  target  extraction  region)  and 
by  an  increase  in  the  RCS  of  the  target. 
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Figure  11.  Bandwidth- 
averaged  RCS 
unprocessed  and 
processed  with  128- 
ramp  ISAR  technique. 
Vertical  lines  show 
RCS  low-noise,  head- 
on  region; 
unprocessed  RCS 
curve  is  higher. 


Azimuth  angle  H 


To  evaluate  the  processing  techniques,  I  calculated  for  each  technique  the 
average  and  the  standard  deviation  of  the  bandwidth-averaged  RCS  in  the 
low-noise,  head-on  region.  The  results  were  then  smoothed  by  a  144-point 
moving-average  filter.  This  filter  was  selected  because  the  192-ramp  ISAR 
image  resulted  in  144  RCS  values  after  clipping.  Therefore,  the  filter  data 
results  should  show  the  trade-offs  between  the  coherent  integration  gain 
achieved  in  ISAR  imagery  and  the  noncoherent  integration  achieved 
through  smoothing.  The  192-ramp  ISAR  image  should  have  less  noise  be¬ 
cause  it  has  coherent  integration  and  a  smaller  extraction  region,  but  it  is 
more  susceptible  to  variations  in  the  measurement  system  over  time.  If 
more  ramps  of  data  were  used,  the  image  would  require  additional  proc¬ 
essing  to  compensate  for  range  migration  and  polar  format  data.  The  re¬ 
sults  are  shown  in  table  3. 

I  calculated  the  RCS  means  by  averaging  (in  meters  squared)  over  aspect 
angle,  then  converting  to  units  of  decibel-square  meter.  The  RCS  standard 
deviations  were  calculated  in  meters  squared,  then  added  and  subtracted 
to  the  mean  v'alue.  The  two  sums  were  converted  to  decibel-square  meter, 
and  then  their  difference  was  calculated  and  divided  by  two.  As  long  as 
the  mean  is  much  larger  than  the  standard  deviation,  the  upper  and  lower 
bounds  calculated  from  the  standard  deviation  are  symmetric  (which  was 
the  case).  This  technique  was  used  in  calculating  the  standard  deviation  in 
the  subsequent  tables. 

The  results  from  table  3  indicate  that  the  ISAR  processing  technique  was 
better  than  the  HRR  processing  technique  and  that  the  number  of  ramps 
used  to  generate  an  ISAR  image  did  not  make  a  significant  difference  in 
the  mean  or  the  filtered  standard  deviation.  The  standard  deviation  was 
slightly  reduced  as  the  number  of  ramps  in  the  ISAR  image  was  increased. 
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The  RCS  results  processed  with  128-ramp  ISAR  images  were  selected  for 
further  analysis.  Figure  12  shows  the  RCS  processed  with  the  128-ramp 
ISAR  technique.  The  data  are  usually  stable,  with  brief  segments  of  high 
variability  (probably  caused  by  the  wind  shaking  the  target). 

For  the  near-broadside  region,  both  the  HRR  and  the  32-ramp  ISAR 
processing  techniques  were  used  in  calculating  the  measured  RCS.  An 
antenna  pattern  correction  was  applied  only  to  the  ISAR  image.  The  RCS 
was  not  averaged  over  aspect  angle,  since  the  SNR  was  much  higher  and 
there  were  significant  variations  in  RCS. 


Table  3.  Standard 
deviation  calculated 
for  low-noise  region 
dominated  by  returns 
from  hemispheres. 


Processing 

technique 

Mean  ±  SD 
(dBsm  ±  dB) 

Smoothed  SD 
(dB) 

ISAR  32 

-14.79  ±0.82 

0.37 

ISAR  64 

-14.85  ±0.67 

0.36 

ISAR  128 

-14.93  ±0.63 

0.36 

ISAR  192 

-14.89  ±  0.58 

0.34 

HRRl 

-14.17  ±  1.85 

0.73 

HRR  2 

-13.45  ±  1.83 

0.94 

Figure  12.  Bandwidth- 
averaged  RCS 
smoothed  with  a  144- 
point  moving-average 
filter. 
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4.  Results 


The  measured  RCS  data  were  processed  and  the  results  compared  with 
simulated  and  theoretical  RCS  predictions.  I  analyzed  the  agreement  for 
two  regions  of  the  target,  but  not  for  a  third,  transitional  region.  As  previ¬ 
ously  discussed,  I  initially  divided  the  target  into  the  broadside  and  the 
head-on  regions.  However,  near  the  transitions  between  these  two  regions, 
both  the  simulated  and  measured  data  had  limitations.  Although  the  RCS 
of  the  measured  data  had  not  increased  significantly,  I  had  to  increase  the 
size  of  the  extraction  region  significantly  to  account  for  nonspecular  scat¬ 
tering.  This  size  increase  lowered  the  SNR.  Further,  the  simulations  could 
not  model  secondary  scattering  effects.  Because  of  these  limitations,  I 
made  no  comparisons  between  measured  and  simulated  data  in  this  tran¬ 
sition  region. 

I  did  analyze  a  limited  broadside  region  and  a  limited  head-on  region.  The 
broadside  region  was  clearly  defined  by  the  RCS  reduction  and  was  lim¬ 
ited  to  approximately  ±3°  from  broadside  of  the  cylinder.  To  avoid  noise 
and  processing  problems,  I  limited  the  head-on  region  to  ±60°  from  the 
apex  of  either  hemisphere;  I  refer  to  this  region  as  the  low-noise,  head-on 
region. 


4.1  Theory  and  Near-Field  Model 

I  performed  a  theorehcal  calculation  of  the  RCS  using  high-frequency,  per¬ 
fect  electrically  conducting,  monostatic  approximations.  The  RCS  (o)  of  a 
sphere  in  the  far  field  is  given  by 

s  =  for  ka>  20 ,  (2) 

where  a  is  the  radius  of  the  sphere  and  k  is  the  propagation  number.  For  a 
4-in.  sphere  at  94.149  GHz  (the  average  usable  radar  frequency),  c  = 
-14.89  dBsm  and  ka  ~  200.  Applying  the  standard  far-field  criterion,  range 
>  2D^/ X,  results  in  a  range  of  26  m.  Since  this  is  smaller  than  the  measure¬ 
ment  range  of  46.6  m,  the  far-field  approximations  should  be  valid. 

At  broadside,  no  theoretical  RCS  could  be  calculated  and  compared  to 
measured  RCS  because  of  several  complicating  factors.  The  cylinder  por¬ 
tion  of  the  target  is  approximately  4.6  m  long  (so  that  this  portion  is  in  the 
near  field),  the  antenna  pattern  is  not  uniform,  and  the  pylons  supporting 
the  target  had  a  significant  effect  on  the  RCS.  However,  I  was  able  to 
develop  a  near-field  model  that  incorporated  these  elements. 

The  near-field  model  for  the  RCS  is  based  upon  the  linear  superposition  of 
electric  fields  induced  by  the  portion  of  the  target  that  was  covered  by  the 
pylons  and  those  induced  by  the  portion  not  covered  by  the  pylons.  The 
model  for  the  pylons  did  not  attenuate  the  electric  field,  but  it  did  produce 
a  phase  shift  in  the  electric  field.  The  modeled  RCS  was  calculated  from 


CT=  E,ii,q)  +  E^ii,q)ei<P 


(3) 
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where  a  =  RCS,  Ef  is  the  electric  field  of  the  target  not  covered  by  the  py¬ 
lon,  Ep  is  the  electric  field  of  the  target  covered  by  the  pylon,  and  (p  is  the 
phase  shift  associated  with  the  pylon.  For  CADDSCAT,  I  applied  equation 
(3)  directly.  For  NcPDT,  I  simulated  the  phase  shift  associated  with  the  py¬ 
lons  by  changing  the  range  to  the  pylons.  For  the  model  results  obtained 
with  CADDSCAT  simulated  data,  I  had  to  apply  an  antenna  pattern  cor¬ 
rection  to  the  measured  data.  On  the  other  hand,  the  model  results  ob¬ 
tained  with  NcPDT  simulated  data  could  be  compared  directly  with  the 
measured  data. 


The  RCS  was  simulated  as  a  function  of  phase  shift.  I  selected  an  appropri¬ 
ate  phase  shift  by  matching  gross  features  of  the  measured  and  modeled 
RCS.  I  theoretically  estimated  the  phase  shift  using 


<P  =  4k 


(~fE-l)d 

A  ' 


(4) 


where  e  is  the  dielectric  constant  of  buoyant  Styrofoam,  d  is  the  thickness 
of  the  Styrofoam,  and  A  is  the  transmitted  wavelength.  According  to  Knott 
et  al  [7],  a  typical  value  for  the  dielectric  constant  of  Styrofoam  is  1.04.  For 
e  =  1.04,  d  =  1/4  in.,  and  A  =  3.2  mm,  the  phase  shift  is  28°.  (The  pylons  also 
have  a  noncoherent  return,  but  it  was  much  smaller  than  the  RCS  of  the 
target,  so  it  was  not  included  in  the  model.) 


4.2  Comparison  of  Simulation,  Theory,  and  Measurement 


4.2.1  Hemispheres 

I  simulated  the  RCS  for  aspect  angles  dominated  by  returns  from  the  hemi¬ 
spheres  using  CADDSCAT  and  Xpatch,  and  compared  the  results  to  theo¬ 
retical  and  measured  values.  The  results  are  presented  in  two  categories: 
unsmoothed  and  smoothed.  Within  these  two  categories  are  both  proc¬ 
essed  measured  data  and  simulated  data.  The  measured  data  consisted  of 
RCS  bandwidth-averaged  statistics  that  were  obtained  with  the  ISAR  proc¬ 
essing  technique  on  128  ramps  of  data.  Simulated  data  consisted  of  RCS 
bandwidth-averaged  data  obtained  with  Xpatch  and  CADDSCAT.  (The 
simulated  data  did  not  have  to  be  processed  for  clutter  and  noise  reduc¬ 
tion.)  Smoothed  data  were  smoothed  with  an  83-point  or  0.5°  azimuth  angle 
moving-average  algorithm. 

Analysis  of  the  data  indicated  that  the  RCS  statistics  that  were  dominated 
by  returns  from  the  hemispheres  should  be  compared  (for  the  simulated 
data)  for  azimuth  angles  between  0°  and  60°,  and  (for  the  measured  data) 
for  the  low-noise  head-on  region  (four  60°  regions).  The  SNRs  for  the 
measured  data  were  higher  at  these  aspect  angles  because  the  antenna 
beam  fully  illuminated  the  target.  Also,  the  location  of  specular  return  did 
not  move  significantly  during  the  formation  of  an  ISAR  image  with  respect 
to  a  range  cell. 

Table  4  shows  the  mean  and  standard  deviation  of  the  bandwidth- 
averaged  RCS.  Figures  13  and  14  show  a  plot  of  the  processed  RCS  values. 
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Table  4.  Measured  and  simulated  RCS  values  for  returns  dominated  by  a  hemisphere. 


Processing  Measured  Xpatch  CADDSCAT  PSM  Theoretical 

(dBsmldB)  (dBsmldB)  (dBsmldB)  (dBsmidB)  (dBsm±dB) 

Unsmoothed  -14.93  ±0.63  -14.931  0.48  -14.90  +  0.14  -14.94  1  0.38  -14.89  1  0 

Smoothed  -14.93  1  0.36  -14.93  1  0.45  -14.90  1  0.13  -14.94  1  0.25  -14.89  1  0 


Figure  13.  Bandwidth- 
averaged  RCS  values 
for  measured  data  and 
a  point  scatter  model. 
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The  amount  of  data  in  the  plots  has  been  reduced  by  a  factor  of  10  to  im¬ 
prove  readability. 

The  results  from  table  4  indicate  that  the  average  RCS  of  each  data  set 
agreed  with  the  theoretical  predicted  RCS  to  within  the  corresponding 
standard  deviation.  The  simulation  run  using  C ADDSCAT  was  able  to  pre¬ 
dict  the  RCS  more  accurately  and  precisely  than  the  simulation  run  using 
Xpatch,  but  significantly  more  computer  processing  time  was  required. 
The  PSM  average  RCS  was  approximately  equal  to  the  average  measured 
RCS,  which  indicates  that  the  processing  had  a  minimal  effect  on  the  sig¬ 
nal.  The  standard  deviation  of  the  PSM  represents  the  processing  noise. 
There  is  a  trade-off  between  the  amount  of  processing  noise  and  the  reduc¬ 
tion  in  measurement  system  noise.  Assuming  independent  white  noise, 
the  signal  to  processing  noise  ratio  for  the  measured  RCS  is  2.4  dB: 

p  0.63^-0.38^ 

SNR  =  10  log,„pil0  log,,, - --2 - =  2.4  dB  ,  (5) 

^  n  U.oo 

where  Pg  is  signal  power  and  is  processing  noise  power. 

4.2.2  Cylinder 

I  simulated  the  RCS  of  returns  dominated  by  the  cylindrical  portion  of  the 
target  using  CADDSCAT  and  NcPDT  for  aspect  angles  near  90°  and  com¬ 
pared  the  results  to  measured  data.  The  measured  data  were  processed 
with  and  without  an  antenna  pattern  correction  included  in  the  calcula¬ 
tions,  as  described  in  section  3.1.4.  The  CADDSCAT  simulation  required  a 
imiform  antenna  pattern  and  NcPDT  used  the  measured  antenna  pattern. 
Correcting  for  the  antenna  pattern  slightly  degraded  the  data,  but  the  ma¬ 
jor  peaks  and  troughs  of  the  data  remain.  I  calculated  the  near-field  RCS  in 
NcPDT  using  the  definition  introduced  by  Wang  [8]: 

Pi  47C 

where  Pj  is  power  transmitted,  P2  is  power  received,  Ci  is  the  near-field 
RCS,  is  the  maximum  antenna  gain,  and  R  is  the  range  to  the  center 
of  the  target.  This  is  equivalent  to  how  the  measured  RCS  values  were 
calculated. 

Initially,  I  processed  the  measured  data  to  determine  the  variation  in  RCS 
over  frequency  for  all  broadside  aspect  angles  except  those  that  were 
highly  attenuated.  These  returns  were  believed  to  be  significantly  affected 
by  the  supporting  pylons,  and  the  RCS  of  these  returns  for  different  fre¬ 
quencies  fluctuated.  1  determined  the  region  for  analysis  by  examining  the 
relative  peaks  in  the  RCS  of  the  measured  data.  Figure  15  shows  the 
bandwidth-averaged  RCS  and  the  RCS  of  a  single  frequency  versus  aspect 
angle.  The  vertical  lines  show  the  aspect  angles  that  were  initially  ana¬ 
lyzed.  The  standard  deviation  of  the  bandwidth-averaged  RCS  was 
7.3  dBsm  and  the  variation  over  frequency  was  5.2  dBsm.  I  calculated  the 
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Figure  15. 

Bandwidth-averaged 
and  single-frequency 
RCS  for  measured 
data. 


latter  by  taking  the  square  root  of  the  average  of  the  RCS  variance  over  fre¬ 
quency  for  each  aspect  angle.  The  result  indicates  that  there  is  more  varia¬ 
tion  over  aspect  angle  than  over  frequency,  suggesting  that  bandwidth 
averaging  is  reasonable.  Another  indication  that  bandwidth  averaging  will 
be  necessary  is  that  the  RCS  model  was  not  detailed  enough  to  exactly  re¬ 
produce  the  constructive  and  destructive  interference  for  each  frequency  at 
each  aspect  angle. 

To  calculate  the  RCS  for  the  average  usable  radar  frequency  of  94.149  GHz, 
I  began  by  using  NcPDT  and  CADDSCAT  to  calculate  the  electric  fields  for 
two  sections  of  the  test  target:  the  section  not  obscured  by  the  pylons  and 
the  sections  obscured  by  the  pylons.  The  RCS  was  then  calculated  by  appli¬ 
cation  of  equation  (3).  A  phase  delay  was  first  inserted  into  the  electric  field 
calculated  for  the  section  obscured  by  the  pylons,  and  the  result  was  then 
superimposed  on  the  electric  field  calculated  for  the  other  section.  Figures 
16  and  17  show  the  effects  of  different  phase  delay  on  the  RCS  simulated 
with  NcPDT  and  CADDSCAT.  The  NcPDT  simulation  included  the  meas¬ 
ured  antenna  pattern  in  the  calculations. 

1  selected  a  reasonable  phase  delay  by  comparing  the  results  from  figure  15 
with  the  results  from  figure  16  and  17: 1  found  the  measured  reductions  in 
the  RCS  at  aspect  angles  that  roughly  corresponded  to  the  locations  of  the 
pylons  and  matched  these  reduced  RCS  values  to  simulated  results.  A 
phase  delay  of  ;r/2  or  90°  was  selected  for  more  detailed  analysis.  This  is 
considerably  larger  than  the  predicted  phase  shift  of  28°,  which  was  too 
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Figure  18.  RCS 
calculated  with 
CADDSCATand 
measured  RCS 
processed  with 
antenna  pattern 
corrections. 


Figure  19.  RCS 
calculated  with 
NcPDT  and  measured 
RCS  processed 
without  antenna 
pattern  corrections. 


small  to  produce  the  variations  found  in  the  measured  data.  (I  could  have 
selected  a  smaller  phase  delay  for  analysis,  but  the  larger  phase  shift  ac¬ 
centuated  the  phenomenon  of  interest.) 

Figures  18  and  19  show  the  measured  bandwidth-averaged  RCS  and  the 
simulated  RCS  for  the  middle  frequency  with  a  phase  delay  of  90°.  In  fig¬ 
ure  18,  an  antenna  pattern  correction  was  applied  to  the  measured  data, 
and  then  a  linear  adjustment  was  performed.  (This  linear  adjustment  was 
required  because  the  original  antenna  pattern  correction  was  not  sufficient. 
Rather  than  readjust  the  antenna  pattern  and  the  data,  I  performed  a  linear 
correction  as  a  function  of  aspect  angle.)  The  adjusted  data  were  used  for 
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further  analysis.  The  results  from  NcPDT  are  also  slightly  off  from  meas¬ 
ured  results,  but  the  difference  is  less.  No  correction  was  applied  to  these 
data. 

Figures  18  and  19  show  that  the  shape  and  the  location  of  the  nulls  in  the 
simulated  data  do  not  correspond  with  those  in  the  measured  data.  The 
lack  of  agreement  on  the  location  of  nulls  indicates  that  the  actual  interac¬ 
tion  with  the  pylons  is  more  complex  than  the  modeled  interaction.  The 
pylons  are  curved  3-D  objects  with  normals  that  are  different  from  those  of 
the  target.  Because  of  this  curvature,  both  the  direction  and  phase  of  elec¬ 
tromagnetic  radiation  propagating  through  the  pylons  could  change  rela¬ 
tive  to  free  space.  These  changes  could  accovmt  for  a  change  in  position  of 
the  RCS  nulls  and  their  shape. 

A  point-for-point  comparison  is  not  appropriate  for  the  data  in  figures  18 
and  19  because  of  the  difference  in  the  nulls  in  the  measured  and  simu¬ 
lated  results,  the  aspect  angle  sensitivity  of  the  data,  and  the  expected  limi¬ 
tations  in  the  target  model.  However,  other  meaningful  quantities  can  be 
analyzed.  Tables  5  and  6  quantitatively  describe  the  agreement  between 
measured  and  simulated  RCS  for  aspect  angles  near  broadside. 

Table  5  shows  the  crossover  angles  at  which  the  major  contributor  to  the 
RCS  chcmged  from  the  cylinder  to  the  hemisphere,  as  well  as  the  mean  and 
standard  deviation  of  the  RCS  between  those  angles.  The  crossover  angles 
were  defined  as  the  angles  at  which  simulated  and  measured  RCS  ex¬ 
ceeded  a  threshold.  The  thresholds  were  chosen  based  upon  the  RCS  of  the 
target  at  broadside  and  effects  of  the  antenna  pattern.  For  the  measured 
data,  the  thresholds  were  4  and  7  dBsm  with  respect  to  crossover  angles  1 
and  2.  For  NcPDT,  the  thresholds  were  3  and  7  dBsm.  For  the  measured 
data  with  an  antenna  pattern  correction  and  CADDSCAT,  both  thresholds 
were  9  dBsm. 

The  measured  and  simulated  results  shown  in  table  5  are  in  agreement. 
The  simulated  and  measured  mean  RCS  and  standard  deviations  are  ap- 


Table  5.  Crossover 
angles  and  RCS 
statistics. 


Data 

type 

Crossover 
angle  1  (°) 

Crossover 
angle  2  (°) 

RCS  mean  ±  SD 
(dBsm±dB) 

Measured 

2.65 

2.69 

10.7  +  2.0 

NcPDT 

2.64 

2.74 

9.7  ±2.6 

Measured  (antenna 

2.70 

2.76 

11.5  ±1.3 

pattern  corrected) 
CADDSCAT 

2.78 

2.72 

11.5  ±0.9 

Table  6.  Mean  and 
standard  deviation  of 
RCS  and  energy 

Data 

type 

RCS  mean  ±  StD 
(dBsm  ±  dB) 

Frequency  mean  ±  SD 
(1/"±1/'’) 

density  spectrum 

Measured 

10.7  ±2.2 

5.6  ±1.3 

statistics  for  data  with 

NcPDT 

9.8  ±1.5 

3.6  ±0.9 

nulls  removed. 

Measured  (no  antenna 
pattern  correction) 

11.7  ±2.0 

5.7  ±1.7 

CADDSCAT 

11.9  ±1.1 

4.0  ±  0.9 
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proximately  equal,  as  are  the  crossover  angles.  The  results  for  the  data 
with  the  measured  antenna  pattern  are  slightly  narrower  than  those  for 
data  without  an  anterma  pattern;  this  difference  is  to  be  expected,  since  a 
narrow  antenna  pattern  would  attenuate  the  RCS  near  the  crossover 
angles.  The  mean  RCS  results  from  NcPDT  are  slightly  lower  and  are  prob¬ 
ably  the  result  of  several  small  errors.  The  antenna  pattern  was  assumed  to 
be  stationary  during  the  simulation,  but  because  the  target  was  offset  from 
the  center  of  the  turntable,  the  anterma  pattern  moved  with  respect  to  the 
target.  This  move  had  a  small  effect  on  the  range  and  aspect  angle.  To  ac- 
coxmt  for  this  effect,  I  adjusted  the  measured  crossover  angles  in  table  5  by 
0.045°. 

Table  6  shows  the  results  of  a  qualitative  analysis  of  the  RCS  and  energy 
density  spectrum  statistics  for  angles  near  broadside  with  the  nulls  re¬ 
moved  from  both  the  measured  and  simulated  data.  The  starting  and  end¬ 
ing  points  for  the  data  were  determined  by  a  threshold  value  that  was  2  dB 
higher  than  the  values  used  to  determine  the  crossover  frequencies  in  table 
5.  The  data  associated  with  a  null  were  characterized  as  starting  from  a 
peak  RCS,  proceeding  to  a  minimum  RCS,  and  ending  at  a  peak  RCS.  The 
mean  and  standard  deviation  of  RCS  and  the  energy  density  spectrum  sta¬ 
tistics  were  calculated  for  these  data.  The  energy  density  spectrum  statis¬ 
tics  characterized  the  mean  and  standard  deviation  of  the  spectrum  of  the 
RCS,  given  by 


(it;)  =  Xzy|S(M;)p Aw  , 


(7) 


B  = 


E(w-(w))^|S(w)|^ 


(8) 


where  I  S(w)  I  ^  is  equal  to  the  energy  density  spectrum.  I  computed  the 
spectrum  by  estimating  the  autocorrelation  function  using  a  triangle  or 
Bartlett  window,  then  performing  a  padded  FFT  [9].  A  noise  threshold  was 
established  in  the  energy  density  spectrum,  and  values  below  the  thresh¬ 
old  were  set  to  zero.  Also,  the  dc  component  in  the  energy  density  spec¬ 
trum  was  set  to  zero.  The  results  are  shown  in  table  6. 


The  results  in  table  6  indicate  that  although  the  model  could  not  reproduce 
the  variations  in  the  RCS  seen  in  the  measured  data,  the  gross  features 
were  simulated.  The  standard  deviation  of  the  RCS  for  the  measured  data 
was  higher  than  that  of  the  modeled  data.  Also,  the  mean  and  standard 
deviation  of  the  energy  density  spectrum  of  the  measured  data  were 
higher  than  for  the  simulated  data.  The  differences  were  not  great,  indicat¬ 
ing  that  the  modeled  results  roughly  corresponded  to  those  of  the 
measurements. 


5.  Conclusion 


The  RCS  of  a  test  target  that  consisted  of  a  cylinder  with  hemispherical 
ends  was  measured  and  modeled,  and  the  results  were  evaluated.  The  tar¬ 
get  was  measured  with  ARL's  W-band,  frequency-agile,  instrumentation 
radar  at  ATC.  The  target  was  placed  on  a  turntable,  and  radar  measure¬ 
ments  were  made  from  -5°  to  365°  azimuth  angle.  The  data  were  cali¬ 
brated,  and  then  a  methodology  was  developed  to  reduce  clutter  and 
noise.  The  radar  return  of  the  test  target  was  simulated  with  CADDSCAT, 
Xpatch,  NcPDT,  and  a  PSM.  These  codes  carmot  predict  scattering  due  to 
traveling  or  creeping  waves,  but  for  most  aspect  angles  these  phenomena 
appear  negligible.  Since  such  secondary  scattering  was  not  significant,  I 
could  eliminate  noise  and  clutter  by  processing  the  measured  data  with 
narrow  target-extraction  regions  that  focused  on  the  specular  return  of  the 
target.  The  methods  used  to  process  and  simulate  the  target  returns  de¬ 
pended  on  the  aspect  angle.  The  target  was  divided  into  two  regions:  head 
on  and  broadside.  The  RCS  of  the  target  in  the  head-on  region  was  domi¬ 
nated  by  radar  returns  from  the  hemispheres,  and  that  of  the  near-broad¬ 
side  region  was  dominated  by  returns  from  the  cylinder.  Subsections  of 
these  regions  were  further  evaluated. 

The  bandwidth-averaged  RCS  statistics  of  the  low-noise  section  of  the 
head-on  region  were  evaluated.  Measured  and  simulated  results  were  all 
in  close  agreement  with  the  theoretical  predicted  RCS.  The  CADDSCAT 
simulation  was  more  precise  than  the  Xpatch  simulation,  but  it  was  also 
more  computationally  intensive.  The  measured  results  demonstrated  that 
these  simulations  could  accurately  extract  RCS  data  from  radar  signals 
with  a  low  SNR.  Specifically,  hemispheres  with  a  theoretical  RCS  of 
-14.9  dBsm  were  measured  in  noise  and  clutter  that  had  an  average  RCS  of 
-3  dBsm.  After  the  measured  data  were  processed,  the  bandwidth- 
averaged  RCS  was  -14.9  dBsm  ±  0.6  dB. 

Analysis  of  the  broadside  region  indicated  that  there  was  an  interaction 
between  the  test  target  and  supporting  pylons.  The  test  setup  was  config¬ 
ured  to  minimize  the  clutter  and  noise  in  the  head-on  region;  as  a  side  ef¬ 
fect,  this  configuration  increased  the  clutter  in  the  near-broadside  region. 
The  target /pylon  interaction  was  simulated  by  an  RCS  model,  for  which 
input  was  generated  by  NcPDT  and  CADDSCAT.  Although  the  model 
could  reproduce  some  of  the  statistics  of  the  measured  data,  there  were  sig¬ 
nificant  discrepancies.  Simulated  and  measured  results  were  in  agreement 
for  the  measured  crossover  angles  of  high-to-low  RCS  and  average  RCS. 
On  the  other  hand,  the  nulls  and  spectral  characteristics  of  the  RCS  be¬ 
tween  the  crossover  angles  were  not  in  close  agreement,  although  they  had 
similar  gross  features.  Because  the  interaction  between  the  target  and  the 
pylons  had  not  been  anticipated,  careful  measurements  were  not  taken 
that  would  allow  detailed  modeling  of  this  phenomenon.  Better  measure¬ 
ments  of  the  setup  might  have  improved  the  model  results.  Future  meas¬ 
urements  should  strive  for  a  pylon  configuration  that  does  not  lead  to  ex¬ 
cessive  interference  at  any  given  aspect  angle,  since  one  never  knows  for 
what  purpose  the  data  will  be  used. 
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